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We report, based on its variation in electronic transport to conpled tension and shear deforma¬ 
tion, a highly sensitive graphene-based strain sensor consisting of an armchair graphene nanoribbon 
(AGNR) between metallic contacts. As the nominal strain at any direction increases from 2.5 to 
10%, the conductance decreases, particularly when the system changes from the electrically neutral 
region. At finite bias voltage, both the raw conductance and the relative proportion of the conduc¬ 
tance depends smoothly on the gate voltage with negligible fluctuations, which is in contrast to that 
of pristine graphene. Specifically, when the nominal strain is 10% and the angle varies from 0° to 
90°, the relative proportion of the conductance changes from 60 to ~90%. 


Graphene has been proposed for many applications due 
to its unique physical properties [IHi , in which the elec¬ 
tronic transport through graphene nanoribbon would be 
affected by line defect Q. Of specific interest to the 
present work, it has been proposed as a strain sensor 
due to ch ang e in the conductivity of graphene-based ma¬ 
terials |6l-[l^ . The field of graphene-based strain sensing 
has rapidly developed since the experimental observation 
of the increase in resistance of CVD graphene samples 
when strain is applied in the direction of the electrical 
current In order to accurately determine the direc¬ 
tion and magnitude of strain, a triaxial graphene-based 
strain sensor composite was proposed; it was found that 
the resistance of graphene may be enhanced or reduced 
by the strain in certain directions 0. The sensitivity 
of graphene to strain originates from the deformation of 
carbon-carbon bonds, which alter the hopping integrals, 
and thus the electronic transport in graphene. 

Though the effect of strain on the band structure of 
graphene and narrow graphene nanoribbons (GNRs) has 
been widely discussed [la - ll^ . strain sensing based on 
electronic transport through graphene and GNRs has re¬ 
cently become of wide interest (Tsl - f^ . In the present 
work, motivated by recent experimental findings of the 
potential benefits of coupled tension/shear deformation 
to simulate the strain generated due to the movement of 
human fingers [3] , we theoretically study coupled tension 
and shear deformation on the transport through armchair 
graphene nanoribbons between metallic contacts. 

Most previous theoretical studies of graphene strain 
sensors have adopted homogeneous junctions. However, 
no lattice mismatch occurs at the interfaces for homoge¬ 
neous junctions, whereas the conductance in unstrained 
GNRs is either zero or one at the Fermi energy E = 0 
[ 2 ^. Unlike most experiments in which the gate voltage 
is only applied to graphene samples, the Fermi energy E 
should vary to investigate electronic transport through 
GNRs. For heterogenous junctions of GNRs between 
quantum wire contacts, transport through GNRs is me¬ 
diated by the gate voltage Vg as in previous experimental 
measurement on electrical properties of graphene samples 
[in. These heterogenous junctions are inspired by the 
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FIG. 1. (Color online.) (a) Schematic illustration of AGNRs, 
connected to two semi-infinite quantum wires. There are N 
and M carbon atoms in x and y directions, respectively. 6 is 
the angle between the direction of applied strain S and a;-axis. 
(b) The tensile and (c) the shear component of the strain at 
30° when the nominal strain is Cna = 10%. 


fact that the contacts are metallic rather than carbon in 
experiments [IHi and the conductance of quantum wire 
contacts is maximal at U = 0 because all channels are 
available to electronic transport. 

However, lattice mismatch may exist at the interfaces 
of heterogeneous junctions. Here we adopt heterojunc¬ 
tions of armchair-edged GNRs (AGNR) between quan- 
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turn wire contacts as discussed in Ref. to minimize 
the effect of lattice mismatch at the interfaces and inves¬ 
tigate the effect of uniaxial plus shear strain on electronic 
transport as illustrated in Fig. [T] The strain is only ap¬ 
plied to the AGNR and impacts the hopping integrals in 
the AGNR. The hamiltonian of the AGNR and contacts 
is described by the tight binding approximation as 

X! “*■ ^3 X! (1) 

where a pair of integers ij indicates the lattice position 
Rij, and Cy (cj^-) is the electron annihilation (creation) 
operator. The summation is over the nearest neighbors 
indicated by (• ■ •). is the hopping integral between 

nearest-neighboring sites indexed by ij and i j . V), is the 
effective gate voltage applied to graphene, which is zero 
in contacts. When Vg slightly varies at the interfaces, the 
conductance does not change much around Vg = Q (^ . 

The deformed configurations of the graphene nanorib¬ 
bons were obtained by molecular mechanics simulations, 
where the strain was obtained via applied displacement 
loading [U, [2^. The rectangular AGNR consisted of 
2832 atoms with a length of L = 10.224 nm and width 
W = 7.018 nm. Displacements were applied in incre¬ 
ments of O.OlA followed by a subsequent energy mini¬ 
mization and relaxation until the change in system en¬ 
ergy was less than 10“^ compared with the previous step. 
The simulations were performed using the open source 
package LAMMPS [2^, and the AIREBO interatomic 
potential with a cutoff of 0.68nm. This potential has 
been shown to accurately describe carbon-carbon inter¬ 
actions resulting in accurate predictions of the mechan¬ 
ical properties of graphene [^. We note that because 
molecular mechanics simulations were performed, which 
are intrinsically at OK, and because all applied displace¬ 
ments were in-plane, there was no out-of-plane buckling 
during the simulation. 

As shown in Fig. [T] (b,c), coupled tension and shear 
were applied onto AGNRs and the corresponding strains 
were calculated as discussed in previous works (2l|. We 
also define the ‘nominal strain’ as the displacement ap¬ 
plied (regardless of the direction) with reference to the 
original length. Once the carbon atomic positions are 
obtained at each value of strain, the hopping along each 
bond (with the length /) [3 

(to = 2.7 eV and a = 0.142 nm) is used as the basis 
for the electronic structure and quantum transport cal¬ 
culations. Due to the 2D nature of our analysis, a bonds 
were not considered in our calculation. 

Strain was applied with a tilted angle 0 (Fig. [T|) from 
0° to 90° at five different angles, i.e. 0°, 30°, 45°, 60° 
and 90°, where 0° loading represents pure tension, 90° 
represents pure shear and the other three couple ten¬ 
sion and shear. To simplify the notation, we also intro¬ 
duce a ‘nominal strain’ Cns, which is defined as the dis¬ 
placement applied over the original nanoribbon length 
regardless of the loading angle. All cases with differ¬ 
ent loading angles are deformed at three stages, namely 



FIG. 2. (Color online.) The current through AGNR depends 
on the bias voltage Vb at Fg = 0 (a) and Vg = to (b). The 
conductance of AGNR changes with Vb at Fg = 0 (lower lines) 
and to (upper symbols) (c), and with the gate voltage in 
AGNR at Vb — 0 and Vb = 0.2to (d). No strain is applied 
(black line) and the nominal strain tns being 10% is applied 
at 0° (red line), 30° (green line), 60° (blue line) and 90° (ma¬ 
genta line) respectively. The size of AGNRs is M = 29 and 
N = 96. Vb and Fg in figures is in the unit of to- 


Cns = 2.5%, 5%, 10%, and we will refer exclusively to 
e„s in the following. The tension {exx) and shear (exy) 
strain components at the three stages for different load¬ 
ing angles are corresponding as: Cxx = 2.5%, 5%, 10% 
and Cxy = 0%, 0%, 0% for 0°; = 2.1%, 4.1%, 8.7% and 

Cxy = 0.5%, 0.9%, 1.8% for 30°; Cxx = 1.6%, 3.3%, 7.8% 
and exy = 0.7%, 1.4%, 2.6% for 45°; exx = 1-1%, 2.4%, 7% 
and exy = 0.9%, 1.8%, 3.1% for 60°; exx = 0%,0%,0% 
and exy = 1.1%, 2.2%, 3.7% for 90°. The deformed 
atomic configuration and the resulting components of the 
tensile and shear strain at 30° when the nominal strain 
is ens = 10% is shown in Fig. [T](b) and (c) [s^. The left 
contact is fixed while the right one is shifted with strain 
similar to Ref. @ when the shear strain is present, and 
the hopping integral in contacts and interfaces are not 
affected by the shear/tension. 

At a finite bias voltage Vb, the current transfer from the 
left contact to the right one is expressed as I{S, Vb,Vg) = 

2®/^ lEf-Vb /2 0’ where S refers to the 

strain and S = 0 stands for no deformation in the 
GNR, e is the electron charge and h is Planck’s constant. 
T{S, E,Vg) is the transmission at the strain S, the en¬ 
ergy E and the gate voltage Vg. Since the effect of the 
gate voltage on the strain sensor has been explored in a 
dual-gate setup we include Vg to estimate the stabil¬ 
ity and applicability of the graphene-based strain sensor. 
Based on the tight binding hamiltonian in Eq.(l) and the 
transfer matrix method, the transmission T(S, E,Vg) is 
obtained through the scattering matrix by solving the 
Schrd ding er equations, which is described in detail in 
Refs. [361. . 
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I{S, Vf), Vg) at 14 = 0 at first increases slowly and then 
sharply with 14 as shown in Fig. [5] (a), since the trans¬ 
mission increases as the GNR deviates from being electri¬ 
cally neutral, which has the lowest density of states [ 2 ^ . 
I{S,Vb,Vg) at Vg = to at first increases linearly and then 
sub-linearly with I 4 as shown in Fig. [2] (b), since the 
transmission decreases as GNR deviates from the high¬ 
est density of states [l^. The conductance is defined as 
G{S,Vb,Vg) = I{S,Vb,Vg)/Vb, and is directly related to 
the transmission at Fermi energy T( S', Ef,Vg) at the limit 
of 14 —>■ 0 [ 13 . The conductance G{S,Vb,Vg) slightly 
increases and decreases at I 4 = 0 and Vg = to respec¬ 
tively when 14 increases, as shown in Fig. me). When I 4 
changes, a large oscillation of the conductance at 14 = 0 
is induced by quantum interference when electrons are re¬ 
flected at the GNR-contact interfaces [H, [ 2 ^ Issl - fs^ , and 
electron-hole asymmetry in conductance is originated by 
odd-numbered rings at the interfaces as shown in Fig. [1] 
[^ . The curve of the conductance versus the gate volt¬ 
age at 14 = 0.2to overlaps the one at 14 = 0 and the 
fluctuation becomes invisible, as shown in Fig. m'i): due 
to the summation of the transmission among the energy 
range [Ef - 14/2, FI/ -|- 14/2]. 

In strained AGNRs, the transport depends on both the 
direction and magnitude of the strain S. When the nom¬ 
inal strength of the strain is 10%, the current at both 
zero and hnite I4 is lower than in undeformed AGNR, 
and gradually increases with the change in angle of the 
applied strain, when the strain varies from 0° to 90° 
combined with shear, as shown in Fig. [2]]a,b). Under 
pure tension at 0°, the chanre in the current completely 
reproduces the observation The slope of the /-I4 
curve (i.e., the conductance) for nearly electrical neu¬ 
tral graphene samples, which corresponds to our result 
at 14 = 0 and which is dependent on the specific ex- 
Mrimental setup, is 1.5 x 10“^U“^ Q and 2 x 10“®O“^ 
0. Our results show that the ballistic conductance of 
the 7.018 nm wide AGNR is around 3 x 10“®O“^ and 
1 X 10“^O“^ at 14 = 0 and Vg = to respectively. The 
conductance dependence on the strain shown in Fig. me) 
is the same as the current dependence on the strain. 

We compare the conductance dependence on pure ten¬ 
sion and pure shear at 90° to estimate the effect of ten¬ 
sion and shear on transport through AGNR. On one 
hand, in contrast to the fact that the conductance at 
—to < Vg < 0.5to under pure tension at 90° is higher 
than undeformed AGNRs as a result of an increase in 
the hopping integrals along the horizontal direction (^ . 
the conductance under pure shear at 90° is always lower 
than that in undeformed AGNRs as \Vg\ > 0.5to except 
for some fluctuations as shown in Fig. [2] Our calcula¬ 
tion is consistent with the observation of the conductance 
dependence on the shear strain Q. On the other hand, 
the conductance around the neutral point (i.e., at small 
jVgl) increases with the angle of tension/shear as shown 
in Fig. me), compared with that under pure tension at 
0°. The maximal conductance of undeformed AGNR be¬ 
ing MGo [Go = 2e^ jh) occurs at E = ±to and I 4 = 0 



FIG. 3. (Color online.) The ratio of the conductance in AG¬ 
NRs before and after application of the strain with the nom¬ 
inal strain being 10% (a) and the nominal strain is 2.5%, 5% 
and 10% at an angle of 0° (b), 45° (c) and 90° (d) between 
a;-axis. The size of AGNR is M = 29 and N — 96. The bias 
voltage is 14 = 0.2to (main panels) and 14 = 0 (insets). 


between graphite contacts and at U = 0 and I 4 = to 
between quantum wire contacts. The maximal conduc¬ 
tance of strained AGNR around I 4 = to decreases as 
a result of deformation in AGNR when tension and/or 
shear is applied to AGNR. The data indicates that elec¬ 
tronic transport through AGNR can be easily mediated 
by the strain when the system deviates from the electri¬ 
cally neutral region. 

In most experiments of graphene-based strain sensors 
0 0 , no gate voltage is applied and graphene may not 
be electrically neutral due to the doping from metallic 
contacts 0 [s^. Recently the effect of the gate voltage 
has been explored and thus we set the gate voltage 
as a variable parameter to provide information such as 
the stability of the strain sensor under different gate volt¬ 
ages. The change of current is usually measured under 
strain and the percentage of the resistance is used to 
estimate the effect of the strain on transport through 
graphene samples 0 0 . Therefore, we use the ratio 
of the conductance compared with that in undeformed 
AGNR, G(S', 14 , 14 )/G(0, 14 , 14 ), to measure the sensi¬ 
tivity of the graphene-based strain sensor as shown in 
Fig. m 

Due to large oscillations in the conductance at 
14 = 0, a large oscillation is also seen in the ratio 
G[S,Vb,Vg)/G{0,Vb,Vg) as shown in the insets of Fig. 
m However, the trends of G[S,Vb,Vg)/G{0,Vb,Vg) are 
still clear. The conductance ratio at 14 = 0.2to shown in 
the main panels of Fig. [3] is relatively smooth at nega¬ 
tive gate voltage, shows a large dip or peak around zero 
gate voltage and slightly decreases as the gate voltage 
becomes more positive. When the nominal strain is 10% 
in Fig. EKa), the conductance slightly increases but is 
smaller than in the undeformed case, as the angle 0 varies 
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from 0° to 90°. As 9 is varied from 0°, 45° and 90° in 
Figs. Eljb-d), the conductance decreases as the nominal 
strain increases from 2.5% to 10%. It is found that the 
conductance shows little change under pure shear at 90° 
as shown in Fig. [H^d). We demonstrate that this kind of 
strain sensor is robust since the relative proportion of the 
conductance is smooth within a wide gate voltage range 


In summary, we have studied a graphene-based 
strain sensor consisting of armchair graphene nanorib¬ 
bon (AGNR) between metallic contacts in response to 
combined tension/shear. The conductance and the rela¬ 
tive proportion of the conductance decreases as the strain 
increases. This kind of strain sensor has relatively higher 
sensitivity to the strength of the strain at finite bias volt¬ 


age and a wide range of the gate voltage when the strain 
is parallel to the armchair edge. 

Finally we comment on the performance of a strain sen¬ 
sor made from a zigzag graphene nanoribbon (ZGNR) 
between quantum wire contacts with a possible lattice 
mismatch at the interfaces. Compared with the case of 
AGNR, the fluctuation of the conductance of ZGNR is 
larger, when the gate voltage changes. The ratio of the 
conductance, G{S,Vb,Vg)/G{0,Vb,Vg), ranges between 
0.8 and 1.4 as \Vg\ < 2to, and the dependence of the 
conductance ratio on the strain are different from that 
seen in Fig. [3] when \Vg\ < to. 
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